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We present an updated analysis of the quark mass dependence of the nucleon mass Mn and nucleon axial-vector 
coupling qa, comparing different formulations of SU(2) Baryon Chiral Effective Field Theory, with and without 
explicit A (1232) degrees of freedom. We discuss the outcome of the corresponding interpolations between lattice 
QCD data and the physical values for these two nucleon observables. It turns out that in order to obtain successful 
interpolating functions at one-loop order, the inclusion of explicit A (1232) degrees of freedom is not decisive for 
the nucleon mass but crucial for qa- A chiral extrapolation of recent lattice results by the LHP collaborations is 
also shown. 


1. THE NUCLEON MASS 

In Ref. ^ we worked out the quark mass depen¬ 
dence of Mjv at the leading and next-to-leading 
one-loop level in manifestly covariant SU{2) 
Baryon Chiral Perturbation Theory (BChPT). 
We neglect isospin breaking effects and trans¬ 
late the functional dependence on the (light) 
quark mass into a pion mass dependence ac¬ 
cording to the Gell-Mann-Oakes-Renner relation. 
We evaluate all the relevant graphs using the so- 
called infrared regularization method [ 2 ] which 
represents a variant of dimensional regularization 
able to treat one-loop integrals involving baryon 
propagators in a way consistent with chiral power 
counting, in a manifestly Lorentz invariant frame¬ 
work. 

In the numerical analysis we have estimated 
parameters not fixed by chiral symmetry by fit¬ 
ting the relevant formulae to a combined set of 
fully dynamical two-flavor lattice QCD data for 
Mn versus obtained by the CP-PACS[21, 
JLQCDlH and QCDSF-UKQCD|5E1 Collabora¬ 
tions. In order to minimize artifacts from dis¬ 
cretization and finite volume effects, we have se¬ 
lected from the set of available simulations those 
with lattice spacing a < 0.15fm and L > 5, 
where L is the spatial size of the lattice. We have 
not considered partial quenching effects, select¬ 


ing lattice calculations with valence quark masses 
equal to sea quark masses. We assume that there 
is an overlap between the region of quark masses 
where SU{2) BChPT is applicable and the range 
of quark masses presently accessible to full-QCD 
lattice calculations, and explore the consequences 
of this assumption. We have restricted ourselves 
to < 600 MeV. 

In Fig^we plot the best-ht curve coming from 
the next-to-leading one-loop expression. We in¬ 
clude the physical point as a constraint. The out¬ 
put parameters show high degree of consistency 
with information from low-energy ttN and NN 
scattering m Statistically compatible results 
have been obtained by fitting with the O(p^) 
Heavy Baryon expression, which coincides with 
the 0{p‘^) expression in the covariant framework 
truncated at in the chiral expansion. This 
confirms the outcome of Ref.|S]. Furthermore, we 
have taken into account finite (spatial) volume 
effects, and fitted to an enlarged set of lattice 
data {L > Ifm) [7] using the 0{p'^) expression 
worked out by the QCDSF-UKQCD collabora¬ 
tion 1^. Our results in the infinite volume are 
very nicely confirmed by those fits 0 Figd also 
shows that the difference between the leading and 
the next-to-leading one-loop result can be kept 
reasonably small over the whole range of m,r that 
we analyze. 
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Figure 1. Solid curve: best fit using the next-to- 
leading one-loop expression for M]^{mTr) 

in covariant SU{2) BChPT, without explicit 
A (1232) degrees of freedom. The physical point 
is included as a constraint. Dashed curve: 
leading-one-loop curve obtained from the 0{p*) 
fit. Dotted curve: 0{p^) result. 


The 0{p'^) pion mass dependence of the pion- 
nucleon sigma term can be obtained by applying 
the Feynman-Hellmann theorem to the expression 
of at the next-to-leading one-loop order. 

According to the outcome of the corresponding 
fit for Mat, we obtain cr^jv = 49 ± 3MeV at the 
physical pion mass pp, in perfect agreement with 
the outcome of the analysis by Gasser, Leutwyler 
and Sainio, = 45 ± 8MeV (H]. 

1.1. Statistical analysis 

In Fig|21we show the statistical error band cor¬ 
responding to the 68% joint confidence region for 
the fit parameters in MAr(m,r) at order |2]. The 
errors on the lattice pion masses have been taken 
into account as well. The inclusion of the physical 
nucleon mass is crucial to shrink the bands below 
about 300 MeV in pion mass. In order to extract 
information from present lattice data about the 
region of small quark masses within our systema¬ 
tic approach, we should either incorporate phe¬ 
nomenological input or perform simultaneous fits 
to several observables characterized by a common 



Figure 2. MAr(m,r): global statistical error band 
associated with the 68% joint confidence region 
for the fit parameters at order p'^. The physical 
point together with the data up to 600 MeV in 
have been included as input. The best-fit curve 
is also plotted. 


subset of low-energy parameters. 

1.2. Including explicit A (1232) degrees of 
freedom 

Following the strategy outlined above we have 
analyzed Mjqim.,:) in the context of an effec¬ 
tive field theory which apart from nucleons and 
pseudo-Goldstone bosons includes the A(I232) 
as an explicit degree of freedom [7|. The theo¬ 
retical framework of our analysis is covariant 
SU{2) BGhPT with infrared regularization in the 
presence of spin-3/2 fields |I()III| . The delta- 
nucleon mass splitting in the SU{2) chiral limit 
is treated as a small parameter and included in 
the power counting according to the so-called 
Small Scale Expansion (SSE)[T2]. Fig|31 shows 
the pion mass dependence of Mat at leading-one- 
loop order in manifestly covariant SSE [7], see 
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Figure 3. MM{mT^) with explicit A(1232) degrees 
of freedom: best-fit curve based on the leading- 
one-loop expression for in manifestly co¬ 

variant Small Scale Expansion. The 68% con¬ 
fidence level error band for the fit at order 
in covariant SU{2) BCHPT with ttN degrees of 
freedom is shown as well. 


Refs.| 10llT| . The curve is again obtained by im¬ 
posing the physical constraint and fitting to the 
available two-flavor fully dynamical lattice QCD 
data compatible with the cuts explained above, 
up to = 600 MeV. As in the previous case, 
the estimated parameters come out of natural size 
and compatible with available information from 
phenomenology |Zj. In FigOlthe best-fit curve is 
plotted together with the 68% error band for the 
0{p^) fit in the scheme with pion and nucleon as 
explicit degrees of freedom. Treating the A(1232) 
as a dynamical variable is not essential for a sati¬ 
sfactory description of the pion mass dependence 
of the nucleon mass. An equally successful inter¬ 
polation can be obtained by “freezing” the delta 
effects into low-energy constants. 


2. THE AXIAL-VECTOR COUPLING 
CONSTANT qa 

We have compared also for gA the chiral effec¬ 
tive field theory schemes with and without ex¬ 
plicit A (1232) degrees of freedom, both in their 
non-relativistic m and manifestly Lorentz in¬ 
variant form with infrared regularization m , see 
also m We have evaluated the contributions 
up to leading and next-to leading one-loop order 
in SU{2) BChPT and used as input for the nu¬ 
merical analysis a set of (quenched) lattice data 
points provided by RIKEN-BNL-Columbia-KEK 
Collaboration m Full-QCD simulations for 
ttItt < 600 MeV are reported by the same collabo¬ 
ration to be fully consistent with the quenched re¬ 
sults In the scheme without explicit A (1232) 
degrees of freedom, both at leading and next-to- 
leading one-loop level, it is impossible to get a 
successful intepolating curve between lattice data 
and the physical point with parameters in agree¬ 
ment with low-energy ttN scattering analyses, cf. 
also the talk by Ulf-G. Meifiner in this work¬ 
shop and Ref. m However, intermediate spin- 
3/2 resonance contributions are known to play 
a very important role in this context. Recall, 
for example, the Adler-Weisberger sum rule m 
which relates the deviation of gA from 1 to pion- 
nucleon dynamics and spontaneous chiral symme¬ 
try breaking: 

2 2fl dq, 

5^ — U-/ — [aT^+p{Lo) — aT^-p{to)]-\-... 

JO ^ 

where the integral of the difference of the 7r+p and 
7r“p total cross-sections is taken over the pion mo¬ 
mentum in the nucleon rest frame and oj denotes 
the pion energy in that frame. Terms of order 
are neglected. The P-wave contributions 
from the nucleon and delta pole terms m give 
total TT^p cross-sections which, according to the 
Adler-Weisberger sum rule, yield gA ~ 1.24. If 
we take into account only the P-wave Born terms 
for static nucleons, we get gA ~ 0.94 d 

We have analyzed gA at the leading one-loop 
order in the framework with explicit A (1232) de¬ 
grees of freedom (non-relativistic SSE) evaluating 
the graphs in Fig 01 The best-fit curve including 
the physical constraint is shown in FigEl The 
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Figure 4. Diagrams contributing to the nucleon 
axial-vector coupling constant gA at leading-one- 
loop order with explicit A (1232) degrees of free¬ 
dom. The wiggly line denotes an external isovec¬ 
tor axial-vector field. 


axial-A-A coupling gi is treated as a free pa¬ 
rameter, while the value of axial-A-A coupling 
is an input for the fit, cf. Fig0 For a choice 
of the latter according to the relativistic expres¬ 
sion of the A ^ ttN strong decay width, the fit 
prefers values of the former close to the SU{4) 
quark model prediction gi = 9/5 g a- We refer to 
EH for an analysis of finite volume effects for gA 
in the same theoretical framework we use. See 
also the contribution by Meinulf Gockeler to this 
workshop. 

In Fig[3 we plot also the leading-one-loop 
(dashed) curve in Heavy Baryon Chiral Pertur¬ 
bation Theory, as extracted from the fit in SSE, 
cf. ■ Clearly, it cannot represent correctly 
the pion mass dependence of gA outside a region 
of extremely small pion masses, as already re¬ 
ported in m The “chiral log” in the leading- 
non-analytic term of the quark mass expansion 
of gA is only visible for pion masses well below 
the physical one m When we analyze the se¬ 
parate contributions to the best fit curve in SSE 
(the solid/dotted one in Fig[2) corresponding to 
the different one-loop diagrams, it turns out that 
the graph with two delta propagators plays an 
important role at the pion masses presently ma¬ 
nageable on the lattice, in order to counterbalance 
the trend shown by the graphs with no delta, see 
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Figure 5. Effects of the inclusion of the explicit 
A(1232) degrees of freedom. The dark solid dots 
refer to simulations by the RBCK Collaboration 

m 


FigEl If the A (1232) is thought to be too heavy 
to propagate, the AA graph starts contributing 
at order in the scheme restricted to pions and 
nucleons only, via a fifth order tadpole with an 
attached external axial field. Therefore the con¬ 
tribution from this diagram starts at higher or¬ 
der with respect to our previous calculation in 
SU{2) BChPT performed up to next-to-leading 
one-loop. 

FigElshows a chiral extrapolation based on the 
leading-one-loop expression in SSE. Here we have 
analyzed very recent lattice data by the LHP col¬ 
laboration m. with two light and one heavier 
flavors, cf. the talk by Wolfram Schroers in this 
workshop. The physical point has not been in¬ 
cluded as input. In the fit an effective low-energy 
coupling has been constrained according to the 
analysis by Fettes et al. of ttN —> ttttN scattering 
, see also m The solid and dashed curves 
correspond to the central, upper and lower value 
for the previously mentioned coupling, as deter¬ 
mined in Ref. m The axial-A-A coupling comes 
out close to the SU{4) quark model prediction. 
Interestingly, the inclusion of a heavier quark in 
the simulations yields results that are statistically 
compatible with the two-flavor case m 
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Figure 6. Solid-curve: best-fit curve at leading- 
one-loop in SSE, FigEl The dot-dashed, dashed 
and dotted curves refer to the diagrams without, 
with one and with two A propagators in Fig^ 
respectively. 
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